ORDER VIA NONLINEARITY IN RANDOMLY CONFINED BOSE GASES 
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A Hartree-Fock mean-field theory of a weakly interacting Bose-gas in a quenched white noise 
disorder potential is presented. A direct continuous transition from the normal gas to a localized 
Bose-glass phase is found which has localized short-lived excitations with a gapless density of states 
and vanishing superfluid density. The critical temperature of this transition is as for an ideal gas 
undergoing Bose-Einstein condensation. Increasing the particle-number density a first-order transi- 
tion from the localized state to a superfluid phase perturbed by disorder is found. At intermediate 
number densities both phases can coexist. 
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1. Introduction 

After the successful realization of the Bose-Hubbard 
model in systems of ultracold atoms in perfectly periodic 
standing wave laser fields [Greiner et ai, 2002] the inter- 
est in the 'dirty boson problem' [Fisher et ai, 1989; Gi- 
amarchi & Schulz, 1987; Giamarchi & Schulz, 1988] has 
also had a strong recent revival [Wang et ai, 2004]. In 
this case cold bosonic atoms in potentials with quenched 
disorder are considered. The disorder appears either nat- 
urally as, e.g., in magnetic wire traps [Folman et al, 
2002; Schumm et at, 2005; Fortagh & Zimmermann, 
2007] or it may be created artificially and controUably 
as, e.g., by the use of laser speckle fields [Dainty, 1975; 
Lye et al., 2005; Clement et al, 2005; Schulte et ai, 
2005] . Historically, the dirty boson problem first arose in 
the context of superfiuid Helium in Vycor [Wong et at, 
1990]. These classical experiments stimulated a number 
of papers in which the Bogoliubov mean-field theory of 
weakly interacting Bose gases was extended to superfluid 
Bose gases moving in a weak disorder-potential [Huang 
& Meng, 1992; Giorgini et al, 1994; Kobayashi & Tsub- 
ota, 2002; Lopatin & Vinokur, 2002; Astrakharchik et 
ai, 2002; Falco et ai, 2007]. In these papers the per- 
turbation and persistence of the superfluid phase in the 
presence of disorder was studied in detail. The prob- 
lem of strong disorder potentials could not be dealt with 
in these theories, however, as this needs more sophisti- 
cated approaches [Hertz et ai, 1979; Navez et al., 2007; 
Yukalov & Graham, 2007; Yukalov et ai, 2007]. Apart 
from this work the contemporary problem of cold atoms 
in quenched random potentials has also been dealt with 
in numerical simulations of lattice systems [Lee & Gunn, 
1990; Ma et al, 1993; Singh & Rokhsar, 1994; Damski 
et al, 2003; Sanpera et al. 2004, Krutitsky et al, 2006]. 
An exception is the one-dimensional case of strongly re- 
pelling 'dirty bosons' which can be solved exactly by a 
mapping, via fermionization, to the problem of the An- 
derson localization of non-interacting fermions [De Mar- 
tino et al, 2005; Krutitsky et ai, 2008]. 



In the present paper we develop an analytical mean-field 
approach to the dirty boson problem of a weakly inter- 
acting Bose gas in a quenched random potential which 
is bounded from below, i.e. we assume that all realiza- 
tions V^(x) are larger than or equal to a fixed lower bound 
Vq. We shall work out our theory for the quite general 
model of a disorder potential V{x.) where the ensemble 
averag e, in the following denoted by • • • , vanishes, i.e. 
V{xi) — 0, and, e.g., the second-order correlation func- 
tion is given according to V{xi)V{x2) — i?^^^(xi, X2). 
Our mean-field theory requires no further assumptions 
about higher correlation functions V{:x.i) ■ . . . ■ V^(x„). In 
the special case of a short-ranged correlation function 
i?^^-'(xi, X2) = i?(5(xi — X2) we find that, in addition 
to the superfluid phase studied in earlier work, for suf- 
ficiently strong disorder there exists a completely local- 
ized phase, either separately, at low number density, or 
in coexistence with a global condensate at intermediate 
number densities. Our approach, from the start, differs 
qualitatively from the earlier mean-field approaches as 
we define and evaluate two basic order parameters of the 
theory from disorder-averages over correlation functions 
of the underlying Bose fields V'*(x, t) and '0(x, r). Here 
T — it stands for the imaginary time which is periodic 
with period hf3, (3 — l/fc^T denoting the reciprocal tem- 
perature. All quantum mechanical expectation values 
here will be understood as being time-ordered with re- 
spect to imaginary time r. For characterizing the super- 
fiuid phase we use the usual condensate density no which 
is defined by the spatial long-range limit of the 2-point 
function according to 



lim (?/'(x, r)-0* (x', T+)) 

|x— x' I — *oo 



no . 



(1) 



Here and in the following the infinitesimally shifted imag- 
inary time r+ = T + 77 with 77 J, is necessary to guaran- 
tee the normal ordering within the underlying functional 
integral representation of the 2-point function. But, in 
addition, we also take into account the possible existence 
of a Bose-glass phase by introducing a corresponding sep- 
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arate order parameter q which is similar to the Edwards- 
Anderson order parameter of a spin glass [Edwards & 
Anderson, 1975]: 



lim |(?/'(x, r)?/;*(x', 

c— x'l — >oc 



'W = {q + no)' 



(2) 



Note that the latter, at T = 0, also follows from the tem- 
poral long-range correlation in direct analogy to the the- 
ory of quantum spin glasses [Read et al., 1995; Sachdev, 
1999]: 



lim ('!/)(x, t)'!/;* (x. 



(3) 



This procedure of introducing a separate Bose-glass or- 
der parameter q is indispensable, we argue here, since 
in a sufficiently large system mutually disconnected local 
mini-condensates of bosons can form at sufficiently low 
temperature. Clearly, if such local, disconnected mini- 
condensates occur, with or without a surrounding sea of 
superfluid bosons, they cannot be described by the usual 
global mean field of a Bose-Einstein condensate, while 
they are well captured by an Edwards- Anderson-like or- 
der parameter. The definition and use of such a Bose- 
glass order parameter in the context of an otherwise very 
familiar model and the discussion of the simple results 
for the Bose-glass phase obtained in this way is the main 
goal of the present work. 

2. Replica Method 

Technically, our approach makes use of the replica 
method which has been established, over the last few 
decades, as a reliable tool for treating disordered sys- 
tems [Mezard et al., 1987; Fischer & Hertz, 1991; Dot- 
senko, 2001]. We start with the functional integral for 
the grand-canonical partition function 



Z 



(4) 



where we employ units with h = 1. The integration is 
performed over all Bose fields ■(/'*(x, r), ■(/'(x, r) which are 
periodic in imaginary time t with period /?. The Eu- 



clidean action in standard notation is given by 



r d 



2M 



+F(x) - ^l] V^x, r) + I r\^. t)V''(x, r)} , (5) 

where M denotes the particle mass, the chemical po- 
tential, and g the strength of the contact interaction. As 
the grand-canonical partition function Z still is a func- 
tional of the disorder potential V{yi), the corresponding 
thermodynamic potential follows from the expectation 
value with respect to disorder 
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■InZ. 



(6) 



In general it is not possible to explicitly evaluate expres- 
sion as the averaging with respect to the disorder 
potential ^(x) and the nonlinear function of the log- 
arithm do no commute: \nZ ^ InZ. An important 
method to perform the averaging procedure prescribed 
by © is provided by calculating the A^th power of the 
grand-canonical partition function Z in the limit A^ — > 0. 
Indeed, from Z^ = 1 + N \nZ + . . . the thermodynamic 
potential ^ is representable as 



1 z^ 

- hm 

/3 N^O 



1 



N 



(7) 



The A^-fold replication of the disordered Bose gas (jH) , ([5]) 
and a subsequent averaging with respect to the disorder 
potential V^(x) can be worked out by using the charac- 
teristic functional of the disorder potential V^(x). Due to 
the above mentioned assumptions about the statistical 
properties of the disorder potential V{x.), its characteris- 
tic functional is of the following form: 



ci • • • 



expji J d^a;j(x)T/(x)| =exp|^^£ J d^Xi 
X / d^x„i?(")(xi,...,x„)j(xi)---j(xn)l . (8) 



I 

The resulting A^-fold replicated imaginary-time action of the weakly interacting Bose gas in a random potential 
reads 



9 1a' 



^„(x,r)-t-| |^„(x,r)in+^i— L (9) 



^W[r,^] ^ f dr f d^xY, |v^:(x,T 

X dri--- dTn I (fxi--- I d^Xn / d^Xn X! "' X! • ■ • ,x„) lV'Qi(xi,ri)l^ • • • l'0Q,^(x„,r„ 

Jo Jo J J J ^^^^ ^^^^ 



where the indices a, a' label the replicated copies of the Bose field. Note that, by definition, all cumulant functions 
i?^"^(xi, . . . , x„) are symmetric with respect to their arguments Xi, . . . , x„. Thus, in leading order n — 2 the random 
potential leads to a residual attractive interaction between the replica fields V'alx, r), iI)ol{'x.,t) which is, in general, 
bilocal in both space and imaginary time. The replicated action Q is not only the starting point for calculating the 
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disorder average of the grand-canonical potential according to ([7]) and 



rN 



(10) 



but also for the calculation of all disorder-averaged correlation functions of ^/'(x, t), 'ip*{'x., t). For instance, the 2-point 
function can be represented within the replica formalism as 



(V'(x,r)V'*(x',T')) = Hm if^ I n <f T^ra' / P^'o' Uo(x, t)C(x', r') e""^'"'!'^* 

^ Q=l [a' = l-' J 



(11) 



3. Hartree-Fock Mean-Field Equations 

Now we shall analyze the action ([5]) in a Hartree-Fock (HF) approximation with respect to the direct and disorder- 
mediated interactions. To this end we split the Bose-fields V'q into a background '^a and fluctuations Sij^a. Due 
to their smallness we only keep terms up to the quartic order in the fluctuations S^Jcn which are then Gaussian 
factorized. With this we obtain the so-called Gross-Pitaevskii equation for the background field, which is modified by 
the disorder-term, in the form 



d 1 

^ - ^ A - ^ -f g E„ (x, r) -f g((5?Aa(x, t) (5^ (x, r+)) \ *„(x, t) 

f-/3 



(12) 



= / dr' [ d^x'f] i?(2)(x,x')|gaa'(x,r;x',r')vI'„'(x',r') + (<5Va'(x',r')<5C'(x',r'+))vI/„(x,T)|+., 
Here we introduced the mean-fields 

Q„„,(x,T;x',t') = «'a(x,T)^';,(x',T') + (5V'a(x,T)5V'a'(x',T')) , Sa(x,T) = (x, T; X, T+ ) , 

where the expectation values are calculated with respect to the effective quadratic action of the fluctuations 
rf3 . N ■ 



d 1 
A - u + 2(7 i;a(x, r) 



/3 

(5?/'q(x, r) — / dr I dr' 



A^'^'^\5ij\5i>]^ dr J d'^xY,^ |'5C(x,t) 
X y d^^y d^x' J2 i?^'Hx,xo|sa(x,T)JC'(x',r')5V'a'(x',r') + Qaa'(x,r;x',r')<5V'a'(x',r')^C(x,T)| 



(13) 



(14) 



Note that the higher-order terms in the ellipsis of 
Eqs. (fT2|) and (fT4|) stem from cumulant functions i?^"^ 
with n > 3. We observe that Eq. (fT5|) defines Qaa' as a 
mean field related to disorder and essentially as the 
total particle number density ria in the replica index a, 
i.e. Eq = ^J^Q. It should be noted that the HF approxi- 
mation leading to and is much simpler but also 
more restrictive than the Hartree-Fock-Bogoliubov ap- 
proximation in which anomalous correlations SipSip and 
dip*S'ip* would also appear. However, not too far below 
the transition temperature the influence of such terms 
is expected to be small [Goldman et al., 1981; Huse & 
Siggia, 1982]. 

In the following we solve our HF approximation for the 
special case of a delta-correlated disorder potential, i.e. 
we assume i?'-^-' (x, x') = i?(5(x— x') from now on. In that 
case the mean-field equations have a replica symmetric 



solutii 



*a(x,T)=v/no , Ea(x, r)=S, 

(3Qa'(x,T;x',T') = Q{t - t') Saa' + q + no (15) 



with the Matsubara decomposition 
1 



(16) 



m— — oo 



and the bosonic Matsubara frequencies Wm = 27rm//3. 
With this ansatz cumulant functions i?^"^ with n > 3 in 
the Eqs. (fT2|) and (fT4| do not contribute in the replica 
limit — > 0. Furthermore, the expectation values in 
(fT2)l and (fTS]) are evaluated as follows. The fluctuation 
action (fT4|) is of the form 
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,S^]^ l^^^drj^^^dr' I d^x I d°x'Y,^ ^ (^^(x, r), <5^„(x, r)) ^ (x, r; x', r') (^Jff^/jij 



N N 



where the Fourier-Matsubara transform of the integral kernel G^^/ (x, t; x', t') reads 



a(k,Lj„ 







a*(k, cjm) 
with the quantities 

a(k, w„j) -iuJm + e(k) - /i + 2^1] - RQm , 



b{k,UJm) 

b*{k,uj^) 



1(17) 



(18) 



(19) 



and the free dispersion e(k) = k^/2M. Note that we have already performed in (|18p the replica limit iV ^ which 
eliminates the Hartree contribution of the disorder. The Fourier-Matsubara transform of the corresponding Green 
function follows from algebraically inverting (jlSp which yields in replica limit ^ 0: 



1 







a*(k, w,„) 





6*(k,a;m) 



/ h{\i,UJrn) 

a{k,UJrn)'' 

\ ^ a*(k,w^)2 J 



As the Green function contains expectation values according to 

A<5V^„(x,r)<5C'(x',r')) 



Gaa' (x, t; X , r 



(5V;(x,r)JV„,(x',T')) ; ' 



(20) 



(21) 



we arrive at the expression 



(5Va(x, r)(5Vj*,(x',T')) = .gi(x,T;x',r')5aa/ 

+g2(x,T;x',T'), (22) 



where the respective Fourier-Matsubara coefficients, in 
the replica limit ^ 0, read 



5i(k,cj„ 
52(k, w„ 



1 



13 [-iiUm + e(k) -^i + 2E- RQ„,] 

R{q + no) 5^fi 
[e(k) - ^ + 2E - i?Qo]2 ■ 



,(23) 
(24) 



Evaluating ([13]) together with p^ " (|24|) determines the 
Matsubara coefficients 



(25) 



-2^ 



3/2 



2E 7ri?2 



and the Bose-glass order parameter 

y/l^R{q + m) (M 



y/-fi 2S - RQo \2tt 



3/2 



(26) 



Inserting these results in (|T2)l and (|13p and taking the 
replica limit ^ we find the Gross-Pitaevskii equation 



(^ - 2i; -I- guQ + RQo) = 



(27) 



and the particle density 



n 

n = - S = g + no + lim V ^e*""" . (28) 
g vio ^ p 

m— — oo 



The remaining Matsubara sum in (|28p is evaluated to- 
gether with (|25p by using the zeta function regularization 
method [Kleinert, 2004]. Applying both the Schwinger 
trick 



1 



1 



dss^-^e-''" 



(29) 



a'' n^) JO 

with the Gamma function F and the Poisson sum formula 



m— — oo 

yields the result 



E 



(30) 



E (-*^-™ + «) = ^.r(-.) 

m— — oo 

with the standard polylogarithmic function 

oo 

u^) = E 



(31) 



n'^ 



(32) 
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FIG. 1: Isotherm in the {fl, gn)-pla.ne with /i = /i — 
■kR^ {M/2-k)^ for two values of the disorder parameter R. 1: 
normal Bose gas; 2a, 2b: superfluid with the dimensionless 
disorder strength y^R{Mp^^^ /2nf^'^ set to 0.2 and 0.48, re- 
spectively; 3: Bose glass. The dimensionless coupling con- 
stant gl{2-K0^!'-^ jMfl'^ is fixed to 1. 



With this we obtain finally the equation of state of the 
system 



1 , 

a 



no 



3/2 



C3/2 ( 



(33) 



with the abbreviation /ir(A*) = iJi"2gn{ii)~'KR^{M /2n)'^ 
for the renormalized chemical potential. Thus, Eqs. (pS)) . 
(P7)) , ([55)) , and pS)) with m = form a closed set of non- 
linear equations which determine the types of order the 
bosons can establish in the underlying random potential 
via their quantum statistics and their interaction. 

4. Phase Diagram 

Let us first consider briefly the superfluid phase defined 
by no ^ 0, as it is well-known from the earlier work 
[Huang & Meng, 1992; Giorgini et ai, 1994; Kobayashi & 
Tsubota, 2002; Lopatin & Vinokur, 2002; Astrakharchik 
et ai, 2002; Falco et at, 2007]. To this end we solve 
(l27)l for the condensate density and use ([25]) for m = 0, 
yielding 



no 



3/2 



(34) 



Then we determine the Bose-glass order parameter from 



^R{M/2tt' 

9\/-fJ'r 



,3/2 



^ _ , M 



3/2' 



(35) 



Inserting all results in ([33)1 we find the equation of state 
for fixed temperature: 



1 



3/2 



3/2' 



H 3 



(36) 



A resulting isotherm is given in Fig. [H where we plot gn 
for two different disorder-strengths as a function of the 



chemical potential p, = fir + '2gn{fir) = M — ■nR?{M /2nY 
both in units of kgT. The isothermal compressibility 
K = n^'^dn/dfi can be read off this graph. The part of 
the isotherm with no 7^ consists of two branches and 
only exists above temperature- and disorder-dependent 
minimal values fl = fim and n ~ rim- The upper branch 
is locally stable, the lower branch is locally unstable. The 
stable branch starts at // = with infinite compressibil- 
ity and, for large number density smoothly approaches 
from below the curve n = [/i - 2,^/ttR{M ^ + 
QtiR^{M /2'kY]/ g. The lower unstable branch likewise 
starts sX fl — jlm with negative infinite compressibility, 
passes through a state of vanishing compressibility and, 
for large jl, approaches from above the line n = fi/2g. 
All this is qualitatively the same as in the case with- 
out disorder, where it is well-known that, in the HF- 
approximation, BEG appears via a first-order transition 
(cf. [Fliefier et ai, 2001] and references therein). The 
latent heat of the transition comes from the additional 
repulsive Fock exchange interaction in the gas and the 
glass phase, which is absent for particles in the conden- 
sate. The new additional aspect here is that the cor- 
responding mean-field HF equations in the presence of 
disorder predict that \tfj\'^ 7^ implies also q 0, i.e. 
that the presence of a global condensate implies also the 
presence of mini-condensates if the potential has a dis- 
ordered component. The inverse is not true, however, as 
we shall now see. 

To this end we turn to the branch of the isotherm with 
no — 0. It consists of the normal Bose-gas branch ob- 
tained by solving ([?7)) with no = 0, then ([^5]) with q = 0, 
and finding the equation of state from ()33p : 



2^) ^3/.(e- 



). 



(37) 



This part of the isotherm ends in a critical point at the 
critical particle number density nc = (M/27r/3c)'^^^C(3/2) 
and chemical potential flc = "^gnc or, equivalently, fir — 
0. In addition, however, we find a Bose-glass part of the 
isotherm (no = 0, q 0) along n = fJ'/2g, i.e. fir = 0, for 
fl > flc, n > nc, which starts in the critical point with 
a continuous phase transition. It turns out to have the 
same critical temperature T° = 2n[n/C{3/2)f^^ /ksM 
as the BEG transition of the ideal Bose gas, where 
Ci'^) = Ci'(l) denotes the Riemann zeta function. This 
means that the interaction has a negligible effect at the 
formation of the independent mini-condensates, which 
make up the Bose-glass phase. Furthermore, it has the 
finite compressibility k — l/2gn^ and is locally stable. 
This part of the isotherm is found by again solving ([27]) 
for no = 0, but now solving ()26p for arbitrary q ^ Q 
by putting ^R{M/2TTf''^ = ^-fi + 2Y, - RQo, which 
corresponds to fj,r — 0, so the equation of state reads 



1 

2^ 



(38) 



and finally solving ([55]) by q 



' nr. The isotherm with 
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the normal part and the Bose-glass part is also indicated 
in Fig.[T] Also indicated there is the value jiM of jl for the 
coexistence between superfluid and normal or Bose-glass, 
respectively, which is fixed by the Maxwell-construction. 
To the left (resp. right) of this value the normal or Bose- 
glass state (resp. superfluid state) are absolutely stable. 
Thus, for the two cases of the disorder strength exhibited 
in Fig. [T] only the case with the larger disorder permits a 
stable Bose-glass state. 

We shall now focus our discussion on properties of the 
Bose-glass phase. As we saw it is characterized by the 
absence of off-diagonal long-range order in space (rto — 0) 
but the presence of long-range order in time (g 7^ 0), 
like in a spin-glass phase. Note that the appearance of 
q spontaneously breaks the continuous J7(l)^-symmetry 
for iV-replicas to a joint C/(l)-symmetry. The latter is 
only spontaneously broken if also no ^ 0. 
To see that the superfluid density vanishes in the Bose- 
glass phase we apply an infinitesimal Galilean boost with 
velocity 5\i, using the Galilei-invariance of the distribu- 
tion of the random potential V{yi), and determine the 
linear response of the momentum density (5g — MrinSu, 
which defines the normal gas density n„. A direct evalu- 
ation yields n„ = (7i(x, r; x, r+) -I- (72(x, r; x, r"*"), thus 
the normal density contains a Bose-glass contribution 
via 52- The superfluid density rig, which is defined by 
Us = n — n„, turns out to be equal to uq which, as we 
have seen, vanishes in the pure localized phase. 
The localization of the Bose-glass states can be inferred 
from the spatial exponential fall-off of the correlation 
function (72(x, r; x', r') describing correlations of the lo- 
cally condensed component. Eq. pi]) allows us to ex- 
tract the temperature-independent localization length 
= 2tt/M^R. Since the HF-approximation is an ef- 
fective free-particle theory, this localization length is in- 
dependent of both the number density n and the inter- 
action strength g. The same localization length is found 
by studying the Anderson localization of a particle of 
mass M in a spatial white-noise potential [Edwards & 
Muthukumar, 1988]. 

From (221) and (P5)) follows the spectral function 
Ai{u;,'k) = -2'<s[gi{Ljjn,'i^)\iuj^=uj-io] of the excitations 



2'kRp{lu) 



2 / \2 



(39) 



with the density of states 



p{uj) 



3/2 



- + Mr 0(W + Mr) 

/TT \ ZTT / 



(40) 



Thus, the excitations are resonances with energies larger 
than -Hr = -p.+2gn+l/2Mil, width 7 = |k|/M^a;, and 
wave numbers larger than the inverse localization length 
S^x- The size —pr of the energy gap in the normal phase 
and the value of the chemical potential fj, = 2gn+l/2MQ 
in the gapless localized phase contain the mean-field shift 
due to the repulsive interaction but also the quantum 




FIG. 2: Schematic phase diagram in the disorder- 
temperature plane where density and s-wave scattering length 
are kept fixed. Thick and thin lines represent first order and 
continuous phase transitions, respectively. 



mechanical localization energy of the particles of mass 
M . The gap —p,r in the normal phase disappears at the 
transition to the Bose-glass phase and remains zero in 
that phase. For w — > the density of states then goes 
to zero as ^/uj. Note that our finding differs qualitatively 
from the lattice case where the density of states remains 
non-zero in the limit w ^ [Fisher et ai, 1989; Ma 
et ai, 1993; Singh & Rokhsar, 1994]. We expect that 
a thorough study of the continuum limit of the lattice 
system would show that this density of states tends to 
zero in that limit and, thereby, permit to confirm our 
result. 

5. Conclusion 

In conclusion we have described a simple mean-field HF 
approach to the weakly interacting disordered Bose-gas. 
Within a single approach this leads to the existence and 
possible coexistence of a superfluid phase and a localized 
Bose-glass phase. To illustrate our results further, let 
us briefly discuss the system also in the temperature- 
disorder (T, _R)-plane at fixed particle-number density 
(see Fig. [2]). At i? = there is a critical temperature 
Tc at which the superfluid (stable for T < Tc) and the 
normal Bose-gas (stable for T > Tc) coexist. Note that, 
due to the weak repulsive interaction with an-^^^ <C 1 
(where a — Mg/Air is the s-wave scattering length), this 
Tc is larger than the critical temperature of the ideal Bose 
gas T^ by about ATc = Tc- T^ « 1.3 an^/^T^ (cf. [Kas- 
tening, 2004] and references therein). Increasing R the 
isotherm in Fig. [1] shifts in such a way (to the right) 
that Tc{R) for fixed n decreases. There is then a critical 
minimal value Rc of R at which — jlc above which 
Tc{R) drops below T°, which is the critical tempera- 
ture for the appearance of the localized Bose-glass phase. 
For stronger disorder the normal high-temperature phase 
upon cooling first makes a continuous transition to the 
localized phase, which is then followed by a first-order 
transition to the Bose-condensed phase. These predic- 
tions should be amply testable in current state-of-the-art 
experiments. 
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